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Recent experiments have emphasized that our understanding of the interplay of electron correlations and randomness in solids 
is still incomplete. We address this important issue and demonstrate that particle-hole (ph) symmetry plays a crucial role in 
determining the effects of disorder on the transport and thermodynamic properties of the half-filled Hubbard Hamiltonian. We 
show that the low-temperature conductivity decreases with increasing disorder when ph-symmetry is preserved, and shows the 
opposite behavior, i.e. conductivity increases with increasing disorder, when ph-symmetry is broken. The Mott insulating gap 
is insensitive to weak disorder when there is ph-symmetry, whereas in its absence the gap diminishes with increasing disorder. 



Introduction: The interplay of disorder and interactions 
is at the heart of many interesting and unexplained phe- 
nomena in condensed matter physics. For example, the 
effects of disorder and interactions in two-dimensional 
(2d) electronic systems acting separately lead to insulat- 
ing behavior of the Anderson and Mott kinds, respec- 
tively. However, experimental findings on silicon metal- 
oxide-semiconductor field-effect transistors (MOSFETs) 
show the occurrence of a conducting phase JlJ, which is 
the result of the combined importance of randomness and 
interactions. Other examples of situations in which 
both disorder and interactions are crucial, yet incom- 
pletely understood, include the formation of local mo- 
ments and the behavior of the susceptibility in doped 
semiconductors B, the superconductor-insulator tran- 
sition and universal conductivity in thin metallic films 
[^|,^), and the pinning of flux lines in type-II supercon- 
ductors. j7| 

In recent years, it has become increasingly clear that 
for non-interacting electrons the presence or absence of 
certain symmetries is crucial in determining the effect of 
disorder on both transport and thermodynamic proper- 
ties, as well as critical properties of the localization tran- 
sition. ||| Recent examples where symmetry considera- 
tions are important are given in the context of quantum 
wires Q and disordered superconductors, Jl(],|ll| where 
chiral, time-reversal, and spin-rotation symmetries play 
an important role. 

In this paper, we examine the effect of different types 
of disorder on both the dynamic and equilibrium ther- 
modynamics of the 2d Hubbard model in the vicinity of 
half- filling, electron density (n) = 1. We demonstrate 
that the presence or absence of particle-hole symmetry 
determines the effect of randomness on the conductivity 
and the Mott gap. 

The Model and Computational Approach: We consider 
the following 2d Hubbard Hamiltonian, 



H = - E *y c U„ - E 44rCk* 



«ik»,<r 



(1) 



Here ty is a bond-dependent hopping matrix element on 
near-neighbor sites (ij), t' ik is a bond-dependent hop- 
ping matrix element on next-near- neighbor sites ((ik)), 
U is an on-site repulsion, and /ij is a site-dependent 
chemical potential. We choose -P(iy) = 1/At for ty € 
[t — A t /2,t + At/2], and zero otherwise, with t = 1 
to set our scale of energy. Similarly, P(i- k ) = 1/AJ 
for t[ k G [t' - A' t /2,t' + A' t /2}, and P(jj, s ) = 1/A M for 
/ij G [— A M /2, +A fl /2], so that the various A measure the 
disorder strength. We will focus on half-filling where the 
effects of interactions are most prominent, as evidenced 
by the formation of antiferromagnetic correlations and a 
Mott-Hubbard charge gap at low temperatures. 

Our computational technique is determinant Quantum 
Monte Carlo (QMC), Jl2| an approach which allows us 
to study much larger numbers of particles than those ex- 
plored with exact diagonalization (and the two-electron 
problem) . |lj];[l4|] In this method the electron-electron in- 
teractions are replaced by a space and imaginary time de- 
pendent Hubbard-Stratonovich field. The integral over 
possible field configurations, which exactly retains the in- 
teractions in the problem, is done stochastically and al- 
lows us to calculate static and dynamic (Matsubara time) 
correlation functions at a fixed temperature T . 

The disordered Hubbard model in Eq. (|l) is particle- 
hole (ph) symmetric when i[ k = /ij =0. That is, under 
the transformation c\ a — > (— l) 1 ^ the Hamiltonian is un- 
changed, and the system is precisely half-filled for all val- 
ues of the parameters in H and also for all T. Therefore, 
while near-neighbor bond and local site disorder both 
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introduce randomness into the system, they differ funda- 
mentally in that site disorder breaks ph-symmetry. 

In order to reveal the effects of disorder quantita- 
tively, we will examine the transport by evaluating the 
temperature- dependent dc-conductivity, 

P 2 



Ode 



-A^(q = 0,T = /3/2) 



(2) 
the 



(with (3 = 1/fcsT) as determined [LJ] from 
current-current correlation function, A xx (q, t) = 
Ox(q,T)i x (-q,0)). Here J x (q,r), the q, r-dependent 
current in the x-direction, is the Fourier transform of, 



c i„c 



From the one- 



3x0) = *E<t *l+a,i(q +x 

electron Green function as a function of imaginary time 
we extract the temperature- dependent density of states at 
the chemical potential N(e = 0) flfj: 



N(0) ~ -/3 G (r = 0, t = /3/2) /vr 



(3) 



These two quantities allow a clear characterization of the 
transport and thermodynamic properties of the system. 

Results: First we discuss the transport properties: in 
Fig. 1, we exhibit the effect of near- neighbor hopping 
(bond) disorder on the conductivity. For all disorder 
strengths A t , at temperatures greater than a character- 
istic temperature related to the Mott gap, the system 
shows metallic behavior with o^c increasing upon lower- 
ing T. The conductivity turns down sharply as the tem- 
perature drops below a characteristic temperature and 
the system shows insulating behavior with txd c decreas- 
ing upon lowering T. In the case of zero randomness, the 
perfect nesting of the Fermi surface in 2d leads to anti- 
ferromagnetic long range order (AFLRO) in the ground 
state with and an associated spin density wave gap for 
arbitrarily small U evolving to a Mott gap at larger U. 
Hopping disorder reduces AFLRO via the formation of 
singlets on bonds with large hopping iy and hence large 
coupling J = tfj/U and ultimately destroys it beyond 
At » l.&t. The fascinating result we have found is 
that insulating behavior in the conductivity nevertheless 
persists to much larger A t . Moreover from the shift of 
the maximum in Fig. 1 we deduce that the mobility gap 
in fact increases with increasing A t . 

The situation is quite different in the case of site dis- 
order, as shown in Fig. 2: at fixed temperature T, as site 
disorder A^ is turned on, the conductivity increases, i.e. 
the Mott insulating state is weakened. jl8| At weak disor- 
der, the conductivity drops with decreasing T, reflecting 
again the presence of the Mott insulating phase. As the 
disorder strength becomes large enough to neglect U , one 
would expect a similar temperature dependence arising 
from Anderson insulating behavior. We believe that in 
all cases the conductivity will ultimately turn over and go 
to zero at low T, but we are limited in these simulations 
to temperatures T > W/48 because of the fermion sign 
problem. Nevertheless, the data for site disorder offer a 
dramatic contrast to that of bond disorder (Fig. 1) where 
randomness decreases the conductivity. 
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FIG. 1. The effect of particle-hole-symmetry preserv- 
ing (near-neighbor) bond disorder in the half-filled Hubbard 
Hamiltonian is to decrease the conductivity Ode- Data is for 
U = 4i on a 8 x 8 square lattice; A f measures the strength of 
the bond disorder. 
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FIG. 2. Canonical site disorder (with strength A M ) en- 
hances the conductivity. Particle-hole symmetric site disor- 
der (with strength A^), as with bond disorder (Fig. 1), sup- 
presses the conductivity. Other parameters are as in Fig. 1. 

What is the underlying reason for the different effects 
of bond and site disorder on the conductivity? There are 
several obvious differences in the effect of bond and site 
disorder on local and even longer range spin and charge 
correlations. Site disorder enhances the amount of dou- 
ble occupancy on the lattice, since the energy cost U of 
double occupancy is compensated by differences in site 
energies. One explanation of why site disorder increases 
Ode is that the concomitant increase in empty sites leads 
to more mobility. This destruction of local moments ul- 



2 



timately also leads to the end of antiferromagnetic order. 
Surprisingly, we find in our simulations that bond dis- 
order has a similar diminishing effect on local moments, 
suggesting that the difference in the behavior of the con- 
ductivity arises from a different origin. 

We argue here that particle-hole symmetry is the uni- 
fying criterion which underlies and determines the effect 
of disorder. As emphasized above, site and bond disor- 
der have rather similar effects on the double occupancy. 
Moreover, the consequences of this effect for a^ c are ex- 
pected to become visible only above a threshold value of 
disorder strength, whereas we observe effects on <7d c al- 
ready for weak disorder. Instead, the key distinction is 
in the presence or absence of ph-symmetry. In order to 
explore this conjecture more fully, we have studied two 
other types of disorder: site disorder that preserves ph- 
symmetry and bond disorder that breaks ph-symmetry 
(by including next-near-neighbor hopping). 

Particle-hole symmetric site disorder is introduced 
by adding random chemical potentials to the Hubbard 
model which couple with opposite sign to the density of 
up and down electrons, i.e. choose ^ = = afXj in 
(|l|) . This type of disorder represents a random (Zeeman) 
magnetic field. For U — ph- symmetric site disorder 
has precisely the same effect as conventional site disorder, 
since moving in a given random chemical potential land- 
scape or one obtained by reversing all the site energies is 
entirely equivalent. However the behavior of the conduc- 
tivity at finite U is dramatically different. Fig. 2 shows 
that ph-symmetric site disorder (with strength AJJ has 
the same effect on Ode as bond disorder, i.e. conductivity 
decreases with increasing A^. 

To seek final confirmation of our conjecture, we have 
also explored the effect of next-near-neighbor (nnn) hop- 
ping and randomness therein. Such longer ranging hy- 
bridization breaks ph-symmetry on a square lattice, since 
it connects sites on the same sublattice. We find that 
such disorder has the same effect as conventional site 
randomness, i.e. increases the conductivity at finite T. 
Thus in all four types of disorder, the behavior of the 
conductivity falls into the appropriate class based on the 
preservation or destruction of ph-symmetry, strengthen- 
ing the evidence that it is this symmetry which is playing 
the crucial role in determining the effect of randomness 
on the transport properties. 

We now turn to thermodynamic properties. The most 
direct measure of the Mott gap is from the compress- 
ibility, or from the behavior of density (n) as a function 
of chemical potential ^t, as shown in Fig. 3. The range 
of fi where {n) is constant (and the system is incom- 
pressible) is a direct measure of the gap in the spectrum. 
Hopping and ph-symmetric site disorder clearly stabilize 
the plateau of the density at half-filling. On the other 
hand, conventional site disorder (with A p = U /2) has 
a compressibility which is indistinguishable (within the 
computational possibilities) from the clean system. 
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FIG. 3. The Mott gap is made more robust by the addi- 
tion of bond disorder or particle-hole symmetric site disorder 
(open and filled squares) of strength A = It — U /2, as indi- 
cated by the response of the density to changes in the chemical 
potential. For canonical site disorder (filled circle) the Mott 
gap is practically unaffected by this strength of randomness. 
[19] Calculations are for T = t/8 = W/64 on a 8 x 8 lattice. 
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FIG. 4. Behavior of the density of states 7r7V(0) at the 
Fermi level and at fixed low temperature as a function of disor- 
der strength A/t for various types of disorder. All data are for 
T — t/6, except data for randomness in next-near-neighbor 
hopping (disorder strength A' t ) which are at temperature 
T = t/5 (the value t' = is used). [19] Other parameters 
are as in Figs. 1 and 2. 

The density of states (DOS) at the Fermi level N(Q) 
gives valuable information on the effect of disorder on the 
Mott gap. In the pure system, QMC studies have shown 
that the DOS exhibits a clear Mott gap with N(0) -> 
as T is lowered to zero. The nonzero values of N(0) 
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we obtain at nonzero T reflect the small residual slopes 
in the plateaus in the (n) vs. fj. plot (cf. Fig. 3); at 
lower T, N(0) approaches zero just as the plateaus be- 
come perfectly flat. The behavior of iV(0) at a fixed low 
T as a function of the strength of the various types of 
disorder is given in Fig. 4. N(0) is rather insensitive 
to ph-symmetric disorder (At and AJJ and is even re- 
duced by it: the Mott gap persists. On the other hand, 
ph-symmetry breaking disorder (A M and A' t ) clearly en- 
hances iV(0), i.e. fills up the Mott gap. 

Our results provide a clear numerical demonstration 
of the key role of particle-hole symmetry. The effects 
can also be understood qualitatively as follows: In the 
clean case, at (n) = 1 and strong coupling, the DOS 
consists of an occupied lower Hubbard band (LHB) and 
an unoccupied upper Hubbard band (UHB), separated 
by a charge gap of the order of U. In the case of ph- 
symmetric disorder, the effect of disorder on LHB and 
UHB is identical. Therefore the Fermi energy remains 
in the middle of the gap: this enables the insulating be- 
havior and Mott gap to stay intact. A stabilized charge 
gap for ph-symmetric site disorder is evident since dou- 
ble occupation is strongly suppressed. For nn-hopping 
disorder a simple argument is less obvious, but the data 
in Fig. 3 clearly show that these two cases fall into the 
same class. When ph-symmetry is broken, the LHB and 
the UHB will be affected differently; different numbers 
of states will appear at either side of the gap. As a con- 
sequence, the Fermi energy ends up in one of the tails 
of the DOS, resulting in an enhanced -/V(O) (cf. Fig. 4) 
and increased conductivity (Fig. 2). The fact that the 
states introduced by disorder are localized [^0| will keep 
the system in an insulating state (cf. Fig. 2). 

Conclusions: In this Letter, we have shown that particle- 
hole (ph) symmetry plays a decisive role in determining 
the effect of randomness on transport and thermody- 
namic properties of the half-filled Hubbard model. By 
exploring four different types of disorder, of which two 
preserve ph symmetry and two break it, we demonstrate 
that a classification by this symmetry allows us to un- 
derstand the effect of disorder on (Tdc(T'), the charge gap 
and the compressibility. The presence of ph symmetry is 
found to have a protective influence on the charge gap. 

A related example where symmetry plays a crucial role 
in the effects of disorder is the case of localization in the 
superconducting phase, where the quasiparticles are de- 
scribed by a Bogoliubov-de Gennes Hamiltonian. In 
this case, one can classify the system according to the 
presence or absence of time reversal and spin rotation 
symmetries, and it is found in one dimension that in the 
absence of spin rotation symmetry, the conductance de- 
cays algebraically with system size, while in the sym- 
metric case it decays exponentially. Therefore, in this 
situation as well, the extra spin rotation symmetry leads 
to a strengthening of insulating behavior. 

The question of the behavior of the half-filled fermion 
Hubbard model as disorder is added is furthermore remi- 



niscent of similar issues in the ph-symmetric boson Hub- 
bard model. H At generic densities, it is believed that a 
new 'Bose glass' phase arises to intervene in the original 
ground state phase diagram between superfluid and Mott 
insulating phases, but the situation at the ph-symmetric 
tip of the Mott lobe is uniquely different. Our work is a 
first step in the analysis of the nature of the behavior of 
the fermionic model. 

The authors thank Piet Brouwer, Andreas Ludwig, and 
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